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ABSTRACT 

We present new spectroscopic and photometric observations of the transiting exo- 
planetary system WASP-3. Spectra obtained during two separate transits exhibit the 
Rossiter-McLaughlin (RM) effect and allow us to estimate the sky-projected angle be- 
tween the planetary orbital axis and the stellar rotation axis, A = 3.3tll de grees. This 
alignment between the axes suggests that WASP-3b has a low orbital inclination rela- 
tive to the equatorial plane of its parent star. During our first night of spectroscopic 
measurements, we observed an unexpected redshift briefly exceeding the expected sum 
of the orbital and RM velocities by 140 m s~^. This anomaly could represent the oc- 
cupation of material erupting from the stellar photosphere, although it is more likely 
to be an artifact caused by moonlight scattered into the spectrograph. 

Subject headings: planetary systems - techniques: radial velocities — stars: individual 
(WASP-3) 
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Introduction 



One of the hottest exoplanets yet discovered is WASP-3b, a giant planet with a mass of 
1.76 Mj that orbits an F7-8V star at a dista nce of only 0.032 AU, giving the planet an equilibrium 
temperature of 1960 K (jPollacco et al.ll2008l ). This close-in, transiting system is an ideal target for 
study with the Rossiter-McLaughlin effect because the host star is both relatively bright {V = 10.5) 
and rapidly rotating (usinij, = 13.4ibl.5 km s~^), and the transit occurs with an impact parameter 
near 0.5. Observations of the Rossiter-McLaughlin effect yield valuable information about the 
system, namely the angle (A) between the sky projections of the planetary orbital axis and the 
stellar rotation axis. This parameter is a basic geometric property of the system and a possible 
clue about the processes of planet formation and orbital mig ration. Detailed descriptions of th e 
Rossiter-McLaughlin effect and its applications can b e found in lOhta et al.l (|2005l ). lGimenea (|2006l ). 
Gaudi &: WinnI (|2007l ). and lFabrvckv &: WinnI (120091) . wh ile ea rly examples of suc h measurements 
for exoplanetary systems are given by iQueloz et al.l (j200d ) and I Winn et al.l (|2005l ) . 



In this paper we report the results of new photometric and spectroscopic observations of WASP- 
3, which were conducted with the primary goal of determining A and the secondary goal of refining 
estimates of other system parameters. In §[5] we describe the observations and data reduction, and 
in §[3] we present the model that was used to fit the data and our results for the system parameters. 
These results are discussed in § |H 



2. Observations and Data Reduction 

2.1. Photometric Measurements 

Photometric observations of transits were conducted on UT 2008 May 15, June 10, June 21, 
and UT 2009 May 12 and May 25, using Keplercam on the 1.2m telescop e at the Fred L. Whipple 



Observatory (FLWO) on Mount Hopkins, Arizona (iHolman et al.l 120061 ). On each night our ob- 
servations spanned the entire transit, except for the first and second nights when we missed the 
transit ingress. For the observations in 2008 we used a Sloan i' filter, and for those in 2009 we 
used a Sloan g' filter. We processed the Keplercam images with standard IRAF procedures for 
bias subtraction and flat-field division. Aperture photometry was then performed on WASP-3 and 
6-20 nearby comparison stars, and the WASP-3 signal was divided by a normalized sum of the 
comparison-star signals. 

We also observed the complete transit of UT 2008 July 13 with the University of Hawaii 2.2m 
(UH 2.2m) telescope on Mauna Kea, Hawaii. We used the Orthogonal Parallel Transfer Imaging 
Camera (OPTIC), which is equippe d with two Lincoln Labs CCID128 orthogonal transfer array 



(OTA) detectors (jTonrv et al.lll997l ). Each OTA detector has 2048x4096 pixels and a scale of 
0'.'135 pixel"^. We took advantage of the charge-shifting capability of the OTAs to create square 
point spread functions (PSF) with side lengths of 39 pixels. This allowed us to collect more light 
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before reaching th e saturation hmit than is possible in normal imaging mode (jHowell et alj 120031 : 



Tonrv et al.ll2005l ). We observed through a Sloan z' filter, and bias subtraction, fi a,t-fielding, and 



apertur e photometry were performed with custom IDL procedures described by iJohnson et al 
(|2009al ). 

Final light curves are shown in Fig. [U after correcting for differential extinction as described 
in § 13.11 The photometric data are given in Table [TJ 



2.2. Radial Velocity Measurements 



We measured the apparent radial velocity (RV) of WASP-3 during the transits of UT 2008 June 19 
and 21 with the High Resolution Echelle Spectrometer (HIRES; Vogt et al. 1994) on the Keck I 
10m telescope at the W. M. Keck Observatory on Mauna Kea. The second of these transits was ob- 
served simultaneously with the FLWO 1.2m telescope as described in the previous section. We also 
observed the system with Keck/HIRES on several other nights in 2008 and 2009, all of which were 
outside of transits except for a single measurement obtained during the transit of UT 2009 June 3. 



Our data collection followed the procedure of I Johnson et al.l (|2008l ) in their measurements of 
HAT-P- 1, and the spectrograph was in the sam e configuration as used for the California Planet 



Search (Howard et al, 



2009 



Marcy et al.l l2005l ). Calibration of the instrumental response and 



wavelength scale was achieved using an I2 gas absorption cell. The slit width was set by the 
0'.'86 B5 decker, and the exposure time ranged from 3 to 5 min, giving a resolution of about 60,000 
at 5500 A and a signal-to-noise ratio of approximately 1 20 pixel^^. Doppler shifts were derived 
from the data using the algorithm of iButler et al.l ()1996l ) with subsequent improvements. For a 
given spectrum, measurement errors were derived from the weighted standard deviation of the 
mean among the solutions for individual 2 A spectral segments. The measurement error ranged 
from 6-10 m s~^. 

The RV data are given in Tableland plotted in Fig.[2j The figures show the expected sinusoidal 
variation outside of transits, due to the line-of-sight component of the Keplerian orbital velocity of 
the star, as well as the shorter-timescale anomalous radial velocity due to the Rossiter-McLaughlin 
(RM) effect that occurs during the transit phase. 

In addition to these two sources of RV variability, an unexpected redshift was observed near 
the egress of the UT 2008 June 19 transit. The RV measured at HJD 2,454,637.001983 is redshifted 
by 140 m ■sT^ with respect to our best-fitting model (described in § 13. 2p . This is much larger than 



both the estimated measurement error of 8.3 m s~ and the scatter of ~15 m s~ between the 
other RV data and the best-fitting model. The previous RV data point, at HJD 2,454,636.991034, 
is not as clearly discrepant but it is the second-largest outlier, with a redshift of 49 m s~^ with 
respect to the best-fitting model. In what follows we will refer to these two consecutive RV data 
points as the "RV spike." They are identified with open symbols in Fig. [2j Casual inspection of 
the spectra constituting the spike revealed nothing unusual. In particular, the visual appearance of 
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Fig. 1. — Photometry of transits of WASP-3, based on observations with the FLWO 1.2m and 
UH 2.2m telescopes. Left. — The data and best-fitting model. Right. — Residuals between the data 
and best-fitting model. 
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the 2-d images, the estimated measurement errors, the signal-to-noise ratio of the spectra, and the 
goodness-of-fit statistic (and other metrics) returned by the Doppler code were all within normal 
ranges. Visual examination of the cross-correlation functions (CCFs) did not reveal anything 
unusual about the spike spectra. The time series of bisector spans was also uninformative, as 
it shows the same pattern that was already observed in the Doppler shifts (the Rossiter effect 
and an additional anomaly during the spike). Nevertheless, we believe the spike is an artifact of 
contamination of the spectra by moonlight, rather than an astrophysical phenomenon, based on 
the following logic: 



1. During our observations, the moon was full and 63° away from WASP-3. Thin cirrus clouds 
were noted at approximately the time of the RV spike, and the mean count rates during the 
"spike" observations declined by approximately a factor of 2. (The exposure times for these 
spectra were increased to compensate for the reduced count rates.) Thus, the spike spectra 
were obtained during an interval of reduced transparency and enhanced sky brightness, leading 
to a possible order-of-magnitude increase in the fractional contamination by moonlight. 

2. The moonlight absorption lines would have appeared at a velocity of approximately +10 km s~^ 
relative to the WASP-3 absorption lines (the difference between the barycentric correction of 
4 km s~^ and the WASP-3 systemic velocity of —5.5 km s~^). This is within the WASP-3 line 
profiles, which are rotationally broadened to 14 km s^^. Given the low contrast (~10%) of 
the WASP-3 lines relative to the continuum, a moonlight contamination of only ~1% of the 
total light would be sufficient to produce a ~100 m s~^ shift in the center-of-gravity of the 
WASP-3 lines. This level of contamination would be visually undetectable in the 2-d images. 

3. T o assess whether 1% moonligh t contamination is reasonable we used the moonlight model 



of iKrisciunas &: Schaeferl (Il99ll ). which takes into account the moon phase, star and moon 
coordinates, and observing bandpass. According to this model, in the absence of clouds the 
moonlight flux would have been ~0.1% of the flux of WASP-3. Since the spike spectra were 
obtained when looking through cirrus, it is plausible that the moonlight fraction increased by 
an order of magnitude to 1%. 



Thus, it seems possible and even likely that a small amount of scattered moonlight was re- 
sponsible for the RV spike. We cannot be certain of this conclusion since we could not devise any 
conclusive statistical test for moonlight at such low levels, given the simultaneous presence of the 
RM effect and the time variability of the instrumental broadening profile and the CCFs. For our 
analysis we omitted the RV spike data from consideration, although we discuss another possible 
interpretation in § HI The investigation of the RV spike was the purpose of the coordinated spec- 
troscopic and photometric observations of the transit of UT 2008 June 21. The RV spike did not 
recur, and no outliers of similar magnitude were observed. 
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Fig. 2. — Apparent radial velocity of WASP-3, based on observations with Keck/HIRES. The 
enlarged (gray) error bars include the "stellar jitter" term of 14.8 m s~^ that was added in quadra- 
ture with the measurement errors. The small (black) error bars represent the measurement error 
estimated internally from each particular spectrum. The open symbols represent the "RV spike." 
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3. Analysis and Results 

3.1. Photometric Model 

The first step in our analysis was to fit a parametric model to the photometric data. We 
simplified the analysis by assuming that the planet's orbital eccentricity is zero, which is consistent 
with previous results (jPollacco et al.l l2008l : iMadhusudhan &: WinnI |2009| ) as well as our own RV 
observations and the expectation that such a close-in orbit has been tidally circularized. The orbit 
was parameterized by the period (P), the midtransit time (T^), the impact parameter (6) in units 
of the stellar radius, th e planet-to-star rad ius ratio (r = Rp/Ri,), and the approximate transit 
duration (T) , defined in ICarter et al.l (j2008l ) as 



T 



P 

Tib 



cos I, 



(1) 



where i is the inclination angle. 



The loss of light during the transit was calculated using the analytic formulae of lMandel &: Ago] 



(j2002l ) for a quadratic limb-darkening (LD) law with coefficients ui and U2- To speed the conver- 
gence of our fitting algorithms, we used the nearly-uncorrelated linear combinations of the LD 
coefficients. 



Un 



ui COS cp + U2 sin c 
U2 cos (p — ui sin ( 



(2) 
(3) 



with (j) = 39° (lPalll2008l ). As a compromise between allowing complete freedom in the LD law and 
placing complete trust in a theoretical LD law, we allowed u'l to b e a free param eter and fixed u'2 
at the value pr e dicted from the PHOENIX atmosphere models of IClaretl (120041 ). Additionally as 



m 



Carter et al.l ((20091), we required ui > and < ui + ti2 < 1, to guarantee that the intensity 



profile of the star decreased monotonically toward the limb. 

In addition, we included two additional free parameters for each night of observations, to 
specify the out-of-transit magnitude (m-oot) and a coefficient (k) for differential airmass extinction. 



mobs = "loot + Am + kz, 



(4) 



where mphs is t h e obs erved magnitude. Am is the calculated magnitude change according to the 
Mandel &: Agoll ()2002l ) formulae, and z is the airmass. As with the LD coefficients, it proved 
advantageous to use fitting parameters that were appropriate linear combinations of moot and k. 

We fit the data from all 6 transits, allowing each of the 6 time series to have independent values 
of Tc, r, k, and moot- We required consistency in b and T, and in the LD coefficients for all those 
data sets observed through the same bandpass. However, we did not require consistency in r, in 
order to check for variations in the transit depth that could be caused by a time- variable pattern of 
starspots or other localized intensity variations across the stellar disk. A time-variable pattern of 
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such irregularities would lead to variations in the transit depth from event to event. For example, 
when there are dark spots located on the visible hemisphere of the star, the fractional loss of light 
during transit is greater than when the star is spotless. These transit-to-transit variations would 
be manifested in our model as an apparent change of the r parameter, even though the actual radii 
of the bodies are not changing. The idea was to seek any evidence for stellar activity, which would 
also provide an alternate explanation for the RV spike. 

To determine the best values of the parameters and their uncert ainties, we used a M arkov 
Chain Moiite Car lo (MCMC) algorithm, implemented as described by iHolman et al.l (|2006l ) and 
Winn et al.l (|2007l ). Specifically we used a Gibbs sampler and tailored our jump sizes so that they 
yielded acceptance rates of approximately 40% for each parameter. The number of links in each 
of our chains was approximately 10^ per parameter, so that the distributions for each parameter 
were well converged. Individual jumps were executed with probability exp(— Ax'i/2), where 



E 



Mobs) 



/j(calc) y 



(5) 



In this expression, Nf is the number of data points, /i(obs) is the observed flux, /j(calc) is the flux 
that is calculated for a given set of model parameters, and cjj is a constant for each light curve, 
determined as follows. First, a value ctq was determined such that when using = do, X/ = -^dof 
for the best-fitting model. (The answer was always close to the standard deviation of the out-of- 
transit data.) Next, we m ultiplied an b y a fa ctor /3 to account for time-correla t ed erro rs, using the 



time-averaging method of iPont et al.l (120061 ) as implemented by IWinn et al.l (l2009al ). Averaging 



times of 10-30 min were used to compute /3, giving results (in chronological order) of 0.92, 1.21, 
1.47, 1.30, 2.45, and 1.67. When /3 was found to be smaller than unity, as in the first case, we used 
(3 = 1, reasoning that time correlations can only increase errors. Finding /3 < 1 represents either a 
statistical fluke or a signal in which the red noise has been underestima. ted by analyzing the model 
residuals rather than the out-of-transit data (see lCarter &: WinnI (|2009l ) for a discussion). 



The MCMC analysis was conducted in several stages, taking advantage of the fact that the 
subset of parameters {Tc, k, moot} are nearly uncorrelated with the other parameters in the model. 
First, we determined the midtransit time (Tc) of each event by allowing Tc, k, and rn-oot to be free 
parameters and fixing all other parameters (which are uncorrelated with Tc) at the values that 
minimize Xf- The results for the midtransit times are given in Tabled along with other previously 
reported midtransit times. 

Our collection of midtransit times was then used to determine the transit ephemeris. The 
results are given in Table [H The linear fit that was used to calculate the ephemeris yielded = 21 
with 7 degrees of freedom, indicating an unacceptable fit. The transit timing residuals from this 
fit are shown in Fig. [3j Given the large value of x^, there are either genuine period variations, or 
the transit time uncertainties have been underestimated in some cases. To be conservative for the 
purpose of planning future observations, we have increased the errors of our ephemeris values by a 
factor of y'21/7 above the formal errors of the linear fit. 
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We determined other system parameters by computing another chain, using our newly deter- 
mined ephemeris, holding Tc, fc, and moot fixed at the best-fitting values for each light curve, and 
allowing all other parameters to vary. The key results are given in Tabled! 

Most notable, as mentioned earlier, is the considerable variation in the planet-to-star radius 
ratio (r), at the 2-3(7 level between the different transits. The mean value of r across all 6 transits 
is 0.1059, from which the individual results differ by as much as 5a. The variation in r can also 
be examined in terms of transit depth (5 = r^). The mean of the 6 5 values is 0.01122, and the 
standard deviation is 0.00080, which is larger than the statistical uncertainty in any individual 
depth measurement. As an estimate of the transit-to-transit variation in the transit depth, we 
sought the value of as such that 



E 



((^i -0.01122) 



5, 



(6) 



with the result as = 0.00075. Thus, the data are consistent with fractional variations in the transit 
depth of order as/5 ^ 7%. 

Unfortunately we cannot be confident in the reality of transit depth variations because of the 
possibility of systematic errors. The light curves shown in Fig. [1] clearly exhibit red noise, especially 
the ^'-band light curves. Our fitting procedure attempts to take the red noise into account but it 
is inevitably imperfect. It is also suggestive that the z'-band light curves gave consistent results for 
(5, while the z' and g' light curves gave larger values for 5; this might be due to systematic errors 
in the treatment of LD. 



Comparison with previous studies reveals th at we found a l arger transit depth than lPollacco et al 



( 20081 ) and, similarly, a larger radius ratio than lCibson et al.l ( 20081 ). Our results for these param- 
eters disagree with each of the two studies by 3(T and So", respectively. Of our other p hotom etric 
parameters, we find agreement with the values tha t were reported by iPoUacco et al.l (j2008l ). but 
less so with those reported by I Gibson et al.l (j2008l ) . The latter authors reported a more precise 



result for h of 0.448 =b 0.014 that disagrees with our result of 0.53lt|};|j|^. They also found a shorter 
transit duration, 2.753lo;o^3 hr, compared to our result of 2.813 it 0.012 hr. The re ason for these 



discre pancies is unclear, but we suspect that at least part of the reason is that the I Gibson et al 
(j2008l ) uncertainties were underestimated because the LD law was held fixed. 



3.2. Radial Velocity Model 

We fitted the RV data with a parametric model that takes into account both the orbital velocity 
of the star (V^) and the anomalous velocity due to the RM effect {/Wr). As in our photometric 
model, we maintained the assumption of a circular orbit. Our orbital model was parameterized by 
the period (P), the midtransit time (Tc), the velocity semi-amplitude {K), and a constant velocity 
offset (7). The offset parameter was needed because the precise Doppler velocities were computed 
with respect to a template spectrum with an arbitrary velocity zero point. Furthermore, as we will 
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Fig. 3. — Transit timing residuals. A linear function of epoch was fit to the transit times given 
in Table O and the calculated times were subtracted from the observed times. Circl es denote mid- 



transit times from this w ork; squares denote midtransit times from earlier studies (jPollacco et al 



2008; 



Gibson et al.ll2008l ). 
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describe shortly, we included an additional velocity offset (A7) specific to the observations that 
were made on UT 2008 June 19. 

The RM effect was parameterized by the projected spin-orbit angle (A) and the line-of-sight 
stellar rotat i onal velocity (wsini^,). To calculate AVr we used the "RM calibration" procedure of 
Winn et al. (j2005l ): we simulated spectra exhibiting the RM effect at different phases of the transit 
and then measured the apparent radial velocity of the simulated spectra^ We found the results to 
be consistent with the formula 



AVr = -{Af)Vp 



1.51 - 0.44 ( ^ 

10.0 km s" 



(7) 



where A/ is the fractional loss of light during the transit and Vp is the line-of-sight component of 
the stellar rotation velocity at the location hidden by the planet. 

Since AVr depends on A/, the RV model is dependent on the photometric parameters. To 
determine these parameters, one could fit the RV and photometric data simultaneously, but we 
found it faster and more convenient to fit the RV data alone using prior constraints on the relevant 
photometric parameters. Specifically we used b, T, and r as parameters in the RV model, the first 
two of which were subjected to Gaussian prior constraints from the photometry. As for r, given the 
discrepancies among the results of fitting individual transits, we adopted a Gaussian prior with a 
central value of 0.1049 and an uncertainty of 0.0037, reflecting the median and standard deviation 
of the 6 different results. 



Thus, our RV fitting statistic was 

2 /ui(obs) -t;i(calc)\^ /6-0.531\^ /r-2.412hr\^ /r- 0.1049\^ 
^^ = Z^[ + — + + nnn.^7 ' 



i=l 



ai J \ ai, J \ ax J \ 0.0037 



where A^'^, is the number of RV data points, Vi{ohs) is the observed RV, t;j(calc) is the RV that is 
calculated for a given set of model parameters, and at and ax are taken from Table H] depending on 
th e sign o f the d ifference in the numerator. The LD coefficients were fixed at the values tabulated 



bv IClaretl (j2004l ) for the PHOENIX model in the g' bandpass (where we find the strongest signal 
from the I2 absorption lines used for RV data calibration). The period and midtransit time were 
fixed at the values determined in § 13.11 



The two measurements constituting the RV spike of UT 2008 June 19, described in § 12.21 
were significant outliers that did not conform to our model, regardless of the parameter values that 
were chosen. Assuming that they were influenced by a physical process that was not described 
by our model — whether moonlight, as we discussed in § 2, or an astrophysical phenomenon — we 



^The starting point of the simulations was a spectrum from AUende-Prietol 1 20021 ) of Procyon (F5, TeS = 6500 K) 
that was chosen because it is a slower rotator than WASP-3 but its spectrum is similar in other respects. For details 
on the construction of the simulated Rossiter- McLaughlin spectra, refer to Winn et al. (2005). 
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omitted these two data points from the fit. We also omitted a third RV point immediately following 
the spike, even though this data point was not an outlier, because it was based on observations 
through cirrus clouds and was thus suspected of being affected by moonlight. (None of the results 
we describe subsequently are altered materially if this third point is instead included in the fit.) 
Furthermore, we allowed the data from UT 2008 June 19 to have an extra velocity offset (A7), in 
addition to the velocity offset (7) applied to all of the RV data. We incorporated this additional 
offset to allow for any systematic offset of the data from that night, again due to moonlight or 
astrophysical phenomenon responsible for the RV spike. 

Taking cjj to be the measurement uncertainty given in Table [21 the minimum Xv '^^^ 1^7 with 
28 degrees of freedom, an unacceptable fit. We enlarged the RV errors by adding 14.8 m s~^ in 
quadrature with the measurement errors in order to achieve a reduced Xv °f unity. The excess scatter 
is referred to as "stellar jitter" and is attributed to intrinsic motions of the stellar photosphere, 
unmodeled orbital motions (due to additional planets or companion stars) and unknown systematic 
errors. This procedure of adding in stellar jitter is common in fitting high-precision RV data, and 
the value of 14.8 m s~^ is in line with previous observations of F stars rotating as rapidly as WASP-3 
(|Saar et al.lll998l ). 



Fig. H] shows the RV data, with the omitted points highlighted, the error bars enlarged to 
encompass the stellar jitter, and the velocity offsets applied. Results from the fit are given in 
Table [H We find the parameters describing the RM effect to be v sin i^, = 14.1+JJ km s"^ and 
A = 3.3I44 degrees. Our re sult for v sin iy, i s in a greement with the value of 13.4 ± 1.5 km s~^ that 
was previously reported by Pollacco et al. (j2008l ). based on an analysis of the spectral line profiles 
rather than on the RM effect. The consistency between these results lends confidence to both 
analyses and our "RM calibration" procedure. Our result for A is new information about this 
planetary system, indicating good alignment between the stellar rotation axis and the projected 
orbit normal. Our result for the velocity semi-amplitude K is 290. Sign m s~ ^, which is about 



3a larger than the value of 251.2ljQ^g m s ^reported by ( Pollacco et al. 2008), but wh ich is in 
agreement with the more recent determination of 274 zt 11 m by ISimpson et al.l ( 20091 ). 



4. Discussion 



We have observed the Rossiter-McLaughlin effect during transits of WASP-3b and found that 
the projected spin-orbit angle is A = 3.31^44 degrees, suggesting that the stellar rotation axis and 
the orbital axis are closely alignedl^ This aligm nent makes WASP-3 siin ilar to the majority of other 
exoplanets who se alignment has been studied (iFabrvckv &: WinnI |2009|) and unlike the misaligned 



systems XO-3 (jHebrard et al.ll2008l : IWinn et al.ll2009bl ). HD 80606 (jPont et al.ll2009l : IWinn et al 



While the first draft of this manuscript was under review, we learned of a similar study of WASP-3 by 

IStg" deg, also suggesting a 



Simpson et"all (|2009l ). Those authors also measured the RM effect, and found A 
relatively close spin-orbit alignment. 
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Fig. 4. — The spectroscopic transit of WASP-3. Top. — The RV data from UT 2008 June 19 
(circles), 2008 June 21 (squares), and 2009 June 3 (triangle). The open symbols show the data 
which were omitted from the fit: the 2 points of the RV spike, and the next data point, which is also 
anomalously redshifted although by a much smaller degree. The constant velocity offsets given in 
Table [5] have been applied and the error bars shown here include stellar jitter. Bottom. — Residuals 
between the data and the best-fitting model. 
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2009cl). WASP-14 (Johnson et alJl2009bl). HAT-P-7 (Narita et alJl2009l : IWinn et alJl2009dl ). Corot 



1 (iPont et alJ 120091 ') ■ and WASP-17 (jAnderson et alJl2oio[ )^ It would seem that for WASP-3, the 



inward migration process did not disrupt the initial coplanarity of the system, or tidal effects have 
coplanarized t he system (although the theoretically calculated rate of coplanarization is too slow 
to be relevant, barker &: Ogilvid ((20091)). 



We also used new RV measurements and transit light curves to determine the system pa- 
rameters, and found agreemen t with previous l y pub lished values, with the exception of the RV 
semi-amplitude K reported bv iPollacco et alj (j2008l ). In addition, we observed an intriguing RV 
anomaly during one of the transits. Our investigation of the spectra led us to the conclusion that 
this "RV spike" was caused by a small amount of moonlight that was admitted to the slit along 
with the light from WASP-3. However, we cannot be completely certain in this interpretation, 
and in any case it is interesting to ask whether there are any astrophysical phenomena that could 
produce such an effect. 

One speculative idea is that the spike represents the planet's passage in front of a starspot. 
When a spot is hidden by a planet, the apparent RV of the star will jump to a larger or smaller 
value, because the RM-like effect of the spot temporarily vanishes. When the spot is no longer 
hidden, the star's apparent RV returns to its expected transiting value. The resulting "spike" in RV 
would be analogous to the photometric "bumps" or "rebr ightenings" t h at hav e bee n seen in transit 



l ight curves and attributed to starspot crossings, as in iRabus et al.l (|2009l ) and iDittmann et al 
(j2009l ). amongst others. 



The time between the RV spike's initial point and its maximum (~ 20 mins) is consistent with 
the ingress or egress duration, as would be required of the spot-occultation hypothesis. However, 
to match the -1-140 m s~^ amplitude of the RV spike, the starspot would need to present a velocity 
contrast with the photosphere, in addition to a possible intensity contrast. This can be understood 
as follows: under the starspot hypothesis the amplitude of the RV spike would have an amplitude 

^v^vA^-^\ (9) 



assuming there is a single spot with radius i?^, intensity J^, and radial velocity Vg that may differ 
from the surrounding photosphere. In the absence of the spot, the material at the spot's location 
would have intensity and radial velocity u^. Because the spike occurred near the transit egress, the 
spot would have been near the receding limb of the star, where = +v sin z^Vl — h"^ ~ +12 km s~^. 
Assuming Rs ~ Rp and setting Aw of Eqn. ([9]) equal to 140 m s~^, we find 

hVs 



-0.06. 



(10) 



Since intensity cannot be negative, the above relation implies that the spot's effective RV must be 
oppositely directed from the local photospheric RV. The planet must cover a blueshifted spot in 
order to produce a redshift as large as 140 m s~^. Specifically, the spot's peculiar RV, defined as 
Vg — Vi,, is approximately —12.7 km s^^ (/^.//s). Unless the spot is very bright {Ig 3> /*) the peculiar 
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velocity just estimated would exceed the sound speed of the photosphere, which we estimate to be 
7 km s~^. The hidden-spot hypothesis for the RV spike thus requires not an ordinary spot, but 
rather an erupting and possibly even supersonic spot. What we are calling a "spot" would actually 
be more akin to an active site undergoing a stellar flare. 

Thus, the RV anomaly could be rationalized as a consequence of vigorous stellar activity, a 
hypothesis which would also produce transit depth variations (see § 3.1). A serious objection to 
this hypothesis, however, is that by the usual measure s WASP-3 is a quiesc ent main-sequence F7-8 
dwarf star for which stellar activity is not expected. iPollacco et al.l (|2008l ) presented no evidence 
for photometric variations, chromospheric emission lines, or any other indication of stellar activity. 
Likewise, there is no detectable Ca II H & K emission in any of our own Keck/HIRES spectra, 
including those obtained during the RV spike. The measured Mt. Wilson S values had an RMS of 
only 0.0013, with no noticeable offset even between the measurements spanning different seasons. 
Considering R'^j^, the ratio of chromospheric emission in Ca II H & K cores to the total stellar 
emission, we find that WASP-3's mean logiJ^j^ = —4.9, giving it an approximate age of 4 Gyr 
(jNoyes et al.lll984l ). By comparison, a very active star such as HD189733 has \o^ E!t^^ = —4.5, 
while a very quiet star such as HD9407 has logi?[jj^ = —4.98, (jWright et al.ll2004l ). Finally, one 
might expect an unusually active star to be an unusually rapid rotator, but the value of v sin i^, = 
14.1 km s~^ is typical of an F7V star. The activity hypothesis would therefore imply a high level 
of activity in a star that one would otherwise expect to be quiescent. 

Although the activity hypothesis is not tenable for WASP-3, we predict that similar observa- 
tions of other systems with more active stars will eventually reveal "RV spikes" when the planet 
occults active regions. Simultaneous photometry and spectroscopy of these events will allow in- 
vestigations of the intensity and velocity contrast of the active regions, relative to the surrounding 
photosphere. Active stars are often excluded from Doppler planet surveys, precisely because the 
excess RV noise caused by stellar activity is a hindrance to planet detection. However, magnitude- 
limited photometric transit surveys have no such selection, and in this regard it is unsurprising that 
planets around ver y active stars are being routinely d iscovered by these surveys, such as CoRoT-2 
(jLanza et al.l l2009l ) and CoRoT-7 (jLeger et al.ll2009l ). The Doppler follow-up for these systems 
will undoubtedly be more difficult than for quiescent stars. But, perhaps these difficulties will be 
compensated to some degree by the prospect that active stars with eclipsing planets will eventually 
enable a deeper understanding of starspots, stellar activity, and stellar flares. 
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Table 1. Relative Photometry of WASP-3 



Heliocentric Julian Date 


Relative Flux 


Uncertainty 


Observatory 


2454601.810140 


0.99857 


0.00134 


1 


2454601.810534 


0.99712 


0.00134 


1 


2454601.810927 


0.99775 


0.00134 


1 


2454601.811321 


0.99772 


0.00134 


1 


2454601.811703 


0.99943 


0.00134 


1 



Note. — (1) Fred L. Whipple Observatory 1.2m telescope, (2) University 
of Hawaii 2.2m telescope. 
(This table is available in its entirety online.) 
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Table 2. Keck/HIRES Doppler Shift Measurements of WASP-3 



Heliocentric Julian Date 


Radial Velocity 
(m s-i) 


Measurement Uncertainty 
(m s-i) 


2454ooo.ooob7i 


Ibz.boolz 


b.4zUoU 


24o4oob.o i o m4 


lUl.z / boz 


b.9zooo 


z4o4Dob.8997oo 


O O CCS'I A A 

5d.D9744 


T "1 r' O CT o 


o /I /I c Qc nnonn^j 
i4o4DoD.yUoyyo 


lob.olzoo 


V 1 QQr\Q 

l.iooVo 


z4o4boD.ylo45z 


lYY.bzZOZ 


6.10846 


2454doo.9oU5o9 


138.14350 


7.13066 


2454bob.94UU9b 


1 1 T CIA OT 

117.ol4o7 


b.o471o 


Z4o4bob.9olzlo 


ol.oo49o 


b.o47oU 


z4o4bob.9oyzoy 


O.zoo 


b.49oUo 


i454DOD.yoyse7 


-39.46792 


7.4UUU2 


2454dod.9oU94z 


-8U. 49891 


6.88382 


2454636.991034 


-63.31213 


7.93881 


2454637.001983 


56.56807 


8.33331 


2454637.014090 


-6.59426 


6.29524 


2454637.060144 


-92.67312 


7.94984 


2454637.090295 


-100.37790 


6.66956 


2454638.816277 


-78.81247 


7.72656 


2454638.825606 


-126.87940 


7.02801 


2454638.836659 


-125.27818 


7.66176 


2454638.845548 


-130.38241 


6.97449 


2454638.854761 


-55.94224 


7.53164 


2454638.865225 


-51.31588 


7.19155 


2454638.877991 


-60.45514 


7.23436 


2454638.890896 


-43.96450 


6.91079 


2454639.006175 


-155.20399 


6.97899 


2454674.740671 


-69.57869 


7.92163 


2454674.744085 


-41.62488 


9.83204 


2454674.868667 


66.23444 


8.56491 


2454674.994730 


182.86504 


9.53239 


2454983.909935 


259.36496 


9.73989 


2454984.890077 


-124.05018 


8.87455 


2454986.041289 


-141.69685 


8.96584 


2454986.920247 


34.36423 


8.97649 



- 21 - 



Table 3. WASP-3 Midtransit Times 



Epoch 


Midtransit time [HJD] 


Reference 


-250 


2,454,143.8503 ± 0.0004 


1 


-2 


2,454,601.86514 ± 0.00027 


2 





2,454,605.55956 ± 0.00035 


3 


12 


2,454,627.72098 ± 0.00031 


2 


18 


2,454,638.80329 ± 0.00031 


2 


30 


2,454,660.96435 ± 0.00021 


2 


59 


2,454,714.52210 ± 0.00036 


3 


194 


2,454,963.84361 ± 0.00072 


2 


201 


2,454,976.77290 ± 0.00051 


2 



Note. References: (1) PoUacco ct al. (2008), 
(2) This work, (3) Gibson et al. (2008). 
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Table 4. WASP-3 Photometric Model Parameters 



Parameter Photometric value 

Orbital Period [d] 
Midtransit time [HJD] 
Impact parameter, b 
Approximate transit duration, T [hr] 
Transit duration [hr] 
Rp/R* (UT 2008 May 15, i') 
Rp/R^ (UT 2008 June 9, i') 
Rp/R^ (UT 2008 June 21, i') 
Rp/R^ (UT 2008 July 13, z') 
Rp/R^ (UT 2009 May 12, g') 
Rp/R^ (UT 2009 May 25, g') 
Limb-darkening coefficient ui {g') 
Limb-darkening coefficient U2 {g') 
Limb-darkening coefficient ui [i') 
Limb-darkening coefficient U2 (i') 
Limb-darkening coefficient ui (z') 
Limb-darkening coefficient U2 {z') 



Note. — The values and uncertainties for P and Tc were deter- 
mined by fitting a straight line to the data given in Table [3] and 
multiplying the formal errors in the fit by -^/x^ = \/3. The tran- 
sit duration is dependent on the radius ratio (Rp/R-^), and hence 
the value quoted is the mean and standard deviation of the 6 tran- 
sit durations determined (using each of the Rp/R-^ values). For 
the other parameters, the quoted value represents the mode of the 
MCMC distribution, and the quoted uncertainties are "la" errors, 
spanning the range between the 15.85% and 84.15% levels of the 
MCMC cumulative distribution. One final note is that the approx- 
imate transit duration (T) refers to the quantity defined in Eqn. [TJ 
making it different from the transit duration. 



1.846834 ± 0.000002 
2454605.55922 ± 0.00019 

n c;oi +0.036 

2 412+0-017 

2.813 ±0.012 

n ims+o-ooo? 
0.1038_o.ooii 
1041+00008 

inn +0 0000 

U.lUii_Q QQ]^^ 

1 099+00006 
1107+0-0014 

U.liU/_Q QQ20 

n int^s+o-oou 
5 96+0 031 

p, 91 c:+0.025 

n 000+0.023 

n 001+0.019 
U.OZi_Q 

Q inQ+0.030 

n OOO+0.022 
U.ZOO_Q 036 
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Table 5. WASP-3 RV Model Parameters 



Parameter Value 



■usini^ [km s"^] 14.lti;3 

A [deg] 3.3lit 

K [m s-i] 290.5l|^ 

Mp [Mj] - 2.041°:°^ 

7 [m s-i] 33.5t^J 

A7 [m s-i] 14.ll6:9 



Assuming = 1.24 M© 



Pollacco et alJ (120081 ') . and 



per 

i = 84.1°, as determined from 
our light curve analysis. 



